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Screening effects from sea pairs on tlie quark-antiquark static potential are analyzed phenomeno- 
logically from the light-quark to the heavy-quark meson spectra. From the high excited light-quark 
meson spectrum a universal form for the screened static potential is proposed. This potential is 
then succesfuUy applied to heavy quarkonia. Our results suggest the assignment of X(4260) to the 
4s state of charmonium and the possible existence of a 5s bottomonium resonance around 10748 
MeV. 



I. INTRODUCTION 



One remaining problem in our understanding of QCD 
has to do with quark confinement in hadrons. We expect 
confinement to be the dominant quark-antiquark {q — q) 
or quark-quark interaction at large separation distances 
and therefore to be determinant to explain the proper- 
ties of highly excited (large sized) hadrons. In recent 
years there has been an important progress in the knowl- 
edge of the spectrum of highly excited mesons in the 
light-quark (u, d) [l[ as well as in the heavy-quark (c, b) 
[2, S| sectors. In particular highly excited light-quark 
mesons show an intriguing hydrogen-like spectral pattern 
Q which can be reproduced, within a non relativistic 
constituent quark model framework, by means of a static 
Q—q interaction which becomes asymptotically coulombic 
[5|. Here we suggest that this asymptotic behavior has 
to do with confinement (non-perturbative giuonic effects) 
and not with the perturbative giuonic coulomb interac- 
tion as suggested in We propose that string break- 
ing gives rise, from a linear confining interaction, to an 
asymptotically coulombic potential. When this screened 
confinement is complemented with a screened "giuonic" 
coulomb interaction an accurate description of the highly 
excited light-quark meson spectrum is achieved. The re- 
sulting static potential depends on four parameters. The 
two of them entering in the confinement term, the string 
tension and the string breaking distance, are considered 
as universal in the sense of having the same values in 
all meson sectors. This is also the case for the orbital- 
angular-momentum parameter related to the onset for 
states of confinement. The remaining parameter, the ef- 
fective "giuonic" coulomb strength, gets different values 
when going from light to heavier quarks. It should be 
pointed out that the potential does not contain any ad- 
ditive constant. The calculated meson masses, obtained 
by adding the mass of the quark and the mass of the an- 
tiquark to the eigenvalues of the Schrodinger equation, 
are directly compared to the experimental meson masses. 
In practice, for equal quark and antiquark masses, the 
quark mass and the effective "giuonic" coulomb strength 
are fixed from two well established experimental meson 
masses in the region of applicability of the static ap- 
proach. 



The successful spectral description obtained in the 
light-quark meson case can be extended to other meson 
sectors, in particular to heavy quarkonia where an ac- 
curate description of the highly excited states may be 
very helpful for an unambiguous quantum numbers as- 
signment. Our results suggest that X(4260) could be the 
4s state of charmonium and that a non-cataloged bot- 
tomonium resonance around 10748 MeV might exist. 

These contents are organized as follows. In Section 
II we establish the general character of the constituent 
quark model approximation we use and the general crite- 
rion of validity of the static potential in the different me- 
son sectors. In Section III the asymptotically coulombic 
potential inferred from the study of the highly excited 
light-quark meson spectrum is derived from a screened 
confinement potential ansatz. The implementation of an 
effective "giuonic" interaction and the consideration of 
an additional L~ dependent correcting factor for it al- 
lows for a precise description of the known static spec- 
trum. In Section IV the same potential is applied to 
heavy quarkonia, charmonium and bottomonium, where 
a distinctive quantum numbers assignment for high ex- 
cited states comes out. Finally in Section V we summarize 
our main results and conclusions. 



II. QUARK MODEL APPROACH 

In the study of the meson spectra, from light to heavy 
quarkonia, we shall rely on a Non Relativistic Con- 
stituent Quark Model (NRCQM) framework. We will 
solve the Schrodinger equation for a static potential. Al- 
though the application of the NRCQM to heavy-quark 
systems, at least for bottomonium, can be taken for 
granted, its application to light-quark systems (r7i„ = 340 
MeV) is always a matter of debate. In the spirit of NR- 
CQM calculations the effective values of the parameters 
take into account, at least to some extent, relativistic 
corrections in the kinetic and potential energies. Actu- 
ally it has been recently shown Q that the known spec- 
trum of highly excited light-quark mesons can be nicely 
reproduced within such a calculation framework despite 
its very relativistic character indicated by the calculated 
values of Pu/rriu > 1. 
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The main distinctive feature of the effective potential 
employed in Q is its asymptotically coulombic tendency. 
Explicitly 

k 

Vlight-quarkir OO) = auTs + Cu (1) 

r 

where cr„ = 932.7 MeV/fm stands for the string ten- 
sion, ku = 2480 MeV.fm for a coulomb strength and 
Cu = 1070 MeV for a constant to fix the origin of the 
potential. The distance represents the onset for the 
screening of the interaction due to the presence of light 
quark-antiquark pairs popping out of the vacuum. The 
value chosen = 1.15 fm is inferred from lattice calcula- 
tions (see and references therein). As for cr„ the value 
used is extracted from the phenomenological analysis of 
the (p, a2, ...) Regge trajectory, see also [6|. Concerning 
the value of fc„ one could tentatively try to ascribe it to 
the chromoelectric one gluon exchange (OGE) interac- 
tion as done in However, as we shall show in the 
next section it may be rather giving account of the long 
distance attenuation, due to string breaking, of the linear 
confining term. 

It should be reminded that the accurate energy de- 
scription of the meson states is linked to the correct long 
distance behavior of their wave functions. Given the rel- 
ativistic character of the fitted spectrum we should not 
trust much, at intermediate and short distances, the non- 
relativistic wave functions obtained from the Schrodinger 
equation. Only very large sized light-quark mesons (in 
our scheme the higher the root mean square (rms) ra- 
dius of the meson the lower the Pq/ruq value), for which 
there are no available data yet, can be considered non- 
relativistic systems. For them the wave function coming 
out from the Schrodinger equation might also be accurate 
at short and intermediate distances. 

For the sake of completeness it is worthwhile to re- 
mind the criterion derived in for the applicability of 
a static potential to a meson {qq) sector within our NR- 
CQM framework. It reads 

< r2 >i/2» 1. (2) 

where < >i/2 stands for the root mean square (rms) 
radius of the meson state. So only for mesons with a large 
size, as compared to 1/mq, the static approach makes 
sense. For u and d quarks (m^ = 340 MeV) this means 

< >^/^>> 0.6 fm. In fact the light-quark meson spec- 
trum has been well reproduced for states with rms radii 
greater or equal than three and a half times this limit: 

< r2 >i/2> 2.1 fm. 

When going to heavier quarkonia we have < 
j,2 >i/2>> 0.4 fm for ss {rus ^ 500 MeV), < >^/^» 
0.14 fm for cc (m^ - 1400 MeV) and < >^/^» 0.04 
fm for bb {nib ~ 4800 MeV) where typical values for the 
constituent quark masses have been chosen. Then, by 
using the same validity factor of 3.5 as in the light-quark 
case, we expect the static approximation to be valid for 



(< r2 >^/^)ss > 1.6 fm, (< r2 >^^^)cc > 0.6 fm and 
(< r2 >i/2)^_ > 0.16 fm. 

Notice though that the constituent quark mass in a me- 
son sector is in our model a parameter to be fixed from 
data and that the values of < >i/2 result from the so- 
lution of the Schrodinger equation. Therefore the estab- 
lished criterion has to be checked a posteriori. Nonethe- 
less its consideration is essential to adequately select the 
specific data to be used to fix the free parameters. Thus, 
for light-quark mesons, data corresponding to L = 4 and 
L ~ 5 states were used since these states are expected to 
have large rms radii due to the presence of the centrifugal 
barrier. 



III. LIGHT-QUARK MESONS 

A. String breaking 

The static q ~q potential can be derived from Lattice 
QCD In the quenched approximation, only valence 
quark qy and antiquark q^ , it has the funnel form 

V(r) =ar~- (3) 
r 

where a is the string tension and C is the strength of 
the coulomb interaction. This potential has to be cor- 
rected at short distances so that C becomes a function of 
r f?\. When including sea quarks, an unquenched poten- 
tial results from the screening of the static sources q^ and 
q^ by light qq pairs created in the hadronic vacuum. A 
parametrization of this effect was proposed twenty years 
ago Unquenched lattice results for the potential be- 
tween two heavy static quarks separated by a distance 
r : ^ 1 f m were described by the potential 

Vscrir)^(<7r-^^ (^-^) 

where /i^^ represented a screening length and C was 
related to the quark-quark-gluon coupling through 

C = (4/3)q;s. The screening factor H{r) = (^^—f^p—^ 

was constructed so that Vscr{r) has a coulombic behav- 
ior at small distances whilst approaching a constant at 
large distances. When applied to heavy quarkonia this 
potential form with effective values of its parameters pro- 
vided a precise description of the spectrum of bb states 
with rms radii smaller than 1.1 fm U[l3|- However, up to 
now, lattice calculations do not allow to extract the pre- 
cise form of the QCD static potential at large distances 
[1, [ll| . So the asymptotic constant behavior should be 
considered as an educated guessing. 

Alternatively an attenuated linear form of confinement 
has been implemented for the asymptotic potential in 
the framework of the QCD string approach (QCDSA) 
that has been successfully applied to light-quark [l3, US] 
and heavy-quark ITsj mesons. In this physical pic- 
ture the string tension a is attenuated at separations 
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r > Ri ~ 1.2 fm becoming a function of r so that for 
r > i?2 = 2.5 fm string breaking occurs with large prob- 
abihty. This attenuation plays a key role to correctly 
obtain the masses of the radial excitations of light-quark 
mesons within this approach. More explicitly the con- 
finement potential reads 

/ X .X A exp{y/^{r - Ri) \ 

VsA{r) = cr{r)r = err 1 - 7— — ^- — 

V B + cxp {^/a[r - Ri) J 

(5) 

with o- = 0.185 GeV^ = 937.5 MeV/fm, 7 ^ 0.4, i?i = 6 
GeV~^ = 1.18 fm and B = 20. The screening factor 
between parenthesis that will be called G{r) henceforth 
varies from ~ 1 for r = to a value of ~ (1 — 7) for 
r > i?2. 

Following the same philosophy as in we shall at- 
tempt to extract information over the variation of the 
confining potential with r from a systematic analysis of 
the meson spectrum within our NRCQM framework. As 
mentioned above the main distinctive feature resulting 
from the application of the NRCQM to the light-quark 
meson spectrum is the coulombic asymptotic behavior of 
the potential as given by Eq. ([T|) . It is then interesting to 
examine the possibility that it may come from confine- 
ment as a result of string breaking. Indeed the form of 
the potential in Eq.((T]) can be derived from the screened 
confinement potential ansatz 

Vconf{r)=<Tr{l-e~'^) = arF{r) (6) 

as can be easily checked by using F{r) ^ — ^ and 
making the identifications (up to order l/r^) 

av = a^rs + C„ (7) 

^ = K (8) 

Then from the numerical values of Um, Cu and fc„ previ- 
ously quoted we get 

V = 2.3 fm (9) 

cr = 925.5 MeV/fm (10) 

Let us realize that the value of a stays within the un- 
certainty interval of the phenomenological string tension 
extracted from the p, 02, ...Regge trajectory as it should. 
Regarding v ~ 2rs note its similarity to R2 — 2i?i in the 
QCDSA. In the same manner u can be interpreted as the 
onset for string breaking to occur with large probability. 
We should realize though that for r > R2 = 2.5 fm G{r) 
keeps an almost constant value whereas F(r) varies in a 
coulombic way. 

It is interesting to establish, from the comparison of 
the spectrum obtained from Vconfii^), Eq.®, with data, 
whether the onset for the states of confinement in our 
model may have been experimentally reached or not. 



These states may be characterized for having a vanishing 
probability of presence for r < r^, being a distance re- 
lated to the confinement scale in QCD. From the quan- 
tum number standpoint this means meson states with 
orbital angular momentum L greater or equal than a 
value Lc- A comparison of the light-quark meson spec- 
trum obtained from VconfiT) with data shows that for 
L — A, for instance, the calculated mass is more than 
100 MeV above the upper limit of the experimental in- 
terval. Although this difference between calculation and 
experiment increases when decreasing L, or equivalently 
decreases when increasing L, the significant discrepancy 
for L = 4 might be suggesting that we are still be far 
from the pure confinement region, i. e. Lc » 4. 

B. Phenomenological static potential 

In order to accurately describe the known meson 
spectra the unquenched confinement potential, Vconfir), 
Eq.®, has to be complemented. The natural way to 
do it is through the incorporation of an effective "glu- 
onic" coulomb interaction so that one recovers at short 
distances, when the effect of qq pairs is negligible, the 
quenched (funnel) form of the potential. For the sake of 
simplicity we shall assume the same screening factor used 
for confinement. Thus the potential reads 

l/,b(r) = (<7r-^)(l-e-^) (11) 
r 

where Xq is the "gluonic" coulomb strength. The 
subindex sb indicates that string breaking has been im- 
plemented in both terms of the potential. Certainly Xq 
keeps some relation with the quark-quark-gluon coupling 
as since the chromoelectric OGE contribution should 
be contained in it. However, corrections to the kinetic 
and potential energies could be also taken into account 
through the effective value of Xq. These corrections may 
include for instance relativistic terms in the kinetic en- 
ergy and in the OGE potential, non-perturbative contri- 
butions to the confinement term and to the quark-quark- 
gluon coupling, etc. Therefore Xq has to be considered 
as a free parameter. To fix it from data we calculate 
the high excited light-quark meson spectrum and require 
that the known states with a high orbital angular mo- 
mentum, L = A for instance, for which we expect the 
static approximation works well, are reproduced. The 
results for Xu — 1065 MeV.fm are presented in Table I. 
We have used for the multiplets the quantum numbers 
notation {L,nr), rir : radial quantum number, as derived 
from the solution of the Schrodinger equation. The mass 
in a multiplet is denoted as M^^n,.- Only states giving rise 
to < >^^'^^ 2.1 fm, for which the static approximation 
makes sense, and for which there are well established ex- 
perimental candidates, are considered. The ordering of 
the states has been chosen to make clear the bias of the 
results: the lower the L the bigger the difference between 
calculated masses and data. 
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(1,4) 
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2219 ±43 
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2.9 


2143 


1947 ± 47 





TABLE L Calculated masses ML,nr and rms-radii < >i/2 
for {L,nr) multiplets from Vsb{r) with m„ = 340 MeV, 
a = 925.5 MeV/fm, = 2.3 fm and A„ = 1065 MeV.fm. 
Experimental average masses as in from references [l[ , 
((M)i,„^)csc , and [l], {{M)L,n^)pDG , are shown for com- 
parison. The superindex f in the (4,1) calculated mass indi- 
cates the average mass value chosen to fix A„. 



This deficiency can be corrected in an ad hoc manner 
by introducing an additional L— dependent factor in the 
"gluonic" coulomb term so that 



Vir) 



ar 



ic(L, + l) ' 



(1 - e"-) iiL^Lc 
(12) 



V{t) = V^onfir) a L>L, 



This form for V{r) satisfies effectively the requirement 
that for L > Lc the contributions to the energy from 
other terms in the potential different than the confining 
interaction Vconfif), Eq.®, are neghgible. Moreover, 
the potential for states with L close below Lc differs lit- 
tle from Vconf{r) as it should. Note also that the lower 
the L the bigger the probability for short quark-antiquark 
separations and the bigger the relativistic corrections to 
the potential and kinetic energies. Therefore the value of 
Xq (corresponding to the strength for L = 0) may be in- 
corporating, at least to some extent, such corrections in 
an effective manner. Actually the same discussion done 
for A„ can be repeated here about the effective charac- 
ter of the "gluonic" coulomb strength A„. Therefore any 
attempt to identify our A„ with the coefficient of the chro- 
moelectric coulomb potential obtained from the OGE in 
QCD is risky. 

The calculated masses from V{r), for A^ = 1600 
MeV.fm and Lc = 16, and their comparison to data, are 
shown in Table II where the states have been now ordered 
according to their sizes. As can be seen the agreement is 
remarkable. We should not forget though that the values 
of the parameters, v, A„ and Lc have been fixed from 
the set of data in Table II. Regarding the other parame- 
ters of the model let us recall that the value of a comes 
from an external input: the phenomenological analysis 
of Regge trajectories involving lowly excited light-quark 
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TABLE IL Calculated masses and rms-radii from V{r) with 
m„ = 340 MeV, a = 925.5 MeV/fm, u = 2.3 fm, A„ = 
1600 MeV.fm and Lc = 16. Notation as in Table I. The 
experimental candidates to be members of the multiplets are 
also indicated. 



mesons. As for r7i„ the value chosen corresponds to the 
average dynamic mass generated by spontaneous symme- 
try breaking in the energy region under consideration, see 
and references therein. 



Let us realize that in most multiplets the difference 
in mass between members of the same multiplet is quite 
small what can be interpreted as the absence of signif- 
icant spin-orbit and tensor contributions for the large 
sized states considered. On the other hand the calcu- 
lated meson states become less relativistic when increas- 
ing {L + Ur). So p/m goes from ~ 1.5 for [L + Ur) = 4 
to ~ 1 for {L + rir) = 6. It is then interesting to give the 
model predictions for higher {L,nr) multiplets for which 
the non-relativistic treatment becomes less effective. The 
average values, from the calculated masses correspond- 
ing to the different {L,nr) combinations giving the same 
{L+fir), are listed in Table III. A look at the Table shows 
the quite small difference with the predictions given in Q 
as could be expected from the same asymptotic behavior 
of the potentials employed and the large meson radii in- 
volved. It should be also added that the limiting mass for 
the light-quark meson spectrum evaluated in [5| remains 
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{L + rir) 

(L, Ur) 


< r' >^l^ 
fm 


MeV 


6 

(1,5), (2,4), (3,3) 
(4,2), (5,1) 


4.6 ±0.8 


2421 ± 11 


7 

(1,6), (2,5), (3,4) 
(4,3), (5,2), (6,1) 


6.4 ±0.7 


2524 ± 8 


8 

(1,7), (2,6), (3,5) 
(4,4), (5,3), (6,2), (7,1) 


8.6 ±0.6 


2592 ± 6 


9 

(1,8), (2,7), (3,6), (4,5) 
(5,4), (6,3), (7,2), (8,1) 


11.0 ±0.5 


2638 ± 5 



TABLE III: Predicted masses and rms radii from V(r) for 
some (L,nr) multiplets with i ± > 6. Parameters as in 
Table II. 



almost unaltered as it is given by 

{Mu^it)uu = m.^ + m-+av = 2809 MeV (13) 

IV. HEAVY QUARKONIA 

The proposed form for the screened potential, V{r), 
should be tested in other meson sectors. Heavy quarko- 
nia, in particular the non-relativistic bottomonium, con- 
stitutes the ideal laboratory to do it since the static ap- 
proximation is expected to be valid for the whole spec- 
trum. To apply V{r) to different meson sectors a crite- 
rion to fix the values of the parameters has to be estab- 
lished. As usual we shall assume that the string tension 
a is quite approximately flavour independent. Regard- 
ing V its value has to do with the screening effect caused 
dominantly by light sea quark-antiquark pairs. Conse- 
quently it seems reasonable to take the same value for 
it in the different meson sectors. The universality can 
be tentatively extended to Lc given its connection to Tc, 
or equivalently to the confinement scale in QCD. Our 
results will justify this extension. As for the "gluonic" 
coulomb strength, \q and the quark mass, m^, they will 
be fixed to get the correct splitting and masses of two 
chosen states. 

Let us remark that we are dealing with a spin inde- 
pendent potential. For s— waves we can assume that the 
experimental energy difference between spin singlet and 
spin triplet states (from V{r) they are degenerate) comes 
mainly from the spin-spin interaction. Taking into ac- 
count that the correction for the spin singlet is in ab- 
solute value three times bigger than for the spin triplet, 
we shall consider our calculated s— wave states to be 
approximately describing spin-triplets. For p and d— 
waves spin-orbit and tensor interactions may give sig- 
nificant contributions to the mass. If we remind that for 



light-quark mesons this kind of contributions were sup- 
pressed for rms radii, as calculated in our model, greater 
than 2 fm, we can expect a similar suppression in heavy 
quarkonia to take place perhaps at shorter distances since 
the strength of the spin-orbit and tensor interactions de- 
creases when increasing the mass of the quark. 



A. Charmonium 

In order to fix Ac and TOc we should rely on the highest 
well established excitations with an unambiguous quan- 
tum numbers assignment. We only have J/'ip{\s), tpi'is) 
and Xc0,cr,c2(lp)- As our model does not contain either 
spin-orbit or tensor interactions that can give account of 
the important mass splitting (140 MeV) in the Ip mul- 
tiplet we should choose J/tp{ls) and V'(2s). However, 
this choice is problematic since the rms radius for J/ip 
obtained from the fixed parameters, < >i/2= 0.4 fm, 
would not satisfy the static condition, (< >^^^)cc 
0.6 fm, showing that the static approximation is invalid 
for J/tp. Instead we shall take for granted the conven- 
tional assignment of ip{AOAO) to V'(3s) and choose ip{2s) 
and tp{3s) as referents to fix the parameters. Notice that 
attributing the mass difference M [ip{2s)] — M [ri{2s)] ~ 
49 MeV to the spin-spin interaction the error in the de- 
termination of M [tp{2s)] due to the non consideration 
of such interaction is only of 12 MeV. For M [tp(3s)] we 
expect even a lower error. 

The results for the static cc spectrum for Ac = 157 
MeV.fm and rric = 1448 MeV are shown in Table IV as 
compared to data through a tentative quantum numbers 
assignment. QCDSA results for the UrL states are also 
shown for comparison. 

As expected the lowest p and d states are not well 
reproduced. For Ip, with < >i/2= 0.7 fm, the dis- 
crepancy goes from 20 MeV for Xc2 to 160 MeV for Xco- 
For Id, with < >i/2= i.o fm, the calculated mass 
differ about 100 MeV from the only known experimen- 
tal candidate. The situation improves extraordinarily for 
< f2 >i/2> 1.5 fm since the calculated masses for the 2d, 
4s, 3d and 5s states can be put in perfect correspondence 
with experimental candidates (the resonances y(4360) 
from Belle and y(4324) from BaBar are assumed to cor- 
respond to the same state). This is a very distinctive 
feature of our model. 

It is noteworthy that X(4260) appears as a natural 
4s state (instead in the QCDSA the 4s state is identi- 
fied with ■0(4415)). Actually the reluctancy to assign 
4s quantum numbers to X{4260) comes to some extent 
from the much higher mass predicted from conventional 
charmonium models fl§\ since experimental data might 
be accommodated by making such a choice [l9| . Let us 
also point out that in our model the F(4660) reported 
only by Belle could also correspond to the overlap of the 
energetically close 7s and 8s states. Our 6s model state 
at 4528 MeV would be missed as well as other ns states 
with n > 9. These excitations would be very close in 
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4039 ± 1 
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4398 
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4335 


4361 ± 18^*^ 
4324 ± 24^" 


4464 
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2.6 


4417 


4421 ± 4 
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4468 




4690 
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4804 
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3.5 


4565 
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4.1 


4611 
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4.4 


4639 


4664 ± 16^" 
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4674 
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0.7 


3574 


Xc2(3556.20±0.09) 
Xci(3510.66±0.07) 
Xco(3414.75 ±0.31) 
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1.2 


3965 


Xc2(3929 ±5±2) 
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4212 






4p 
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TABLE IV: Calculated cc masses and rms radii from V{r). 
The superindex f indicates the masses used to fix Ac — 157 
MeV.fm and rric = 1448 MeV. Masses for experimental candi- 
dates, {ML,n,.)pDG, have been taken from 2l unless otherwise 
stated by means of a superindex: Be for ]3elle data [3, Ba 
for BaBar data JJj]- For p waves we quote the npo, npi and 
np2 states. Masses calculated in the QCD string approach 
[14| . {ML,nr)sA, are also shown for comparison. 

energy what could make difficult their experimental dis- 
entanglement despite the fact that the limiting mass of 
the spectrum is still quite far above 

{Murmt)c-c =i m, + + cjv = 5025 MeV (14) 

It should be also remarked the quite relativistic charac- 
ter of the fitted spectrum with {pc/mcY ^ 0.25 — 0.16. 
Nonetheless the values of the wave functions at the origin 
for ns states {n : 2, 3, 4, 5) differ at most a 15% from the 
ones obtained from the solution of the Salpeter equation 
in the QCDSA Hence quite similar results (within a 
20% of difference) would be obtained for the di-electron 
widths and the same conclusion inferred: the measured 
values for ns states {nd [n : 1,2)) are systematically 
smaller (much bigger) than the calculated ones. This 
can be explained by the presence of s — d states mixing 
as a consequence of their coupling to open channels. This 
mixing would significantly modify the values of the s and 
d wave functions at the origin. On the other hand the 
very good fit obtained for the spectrum without mixing 
suggests that this should not have any significant effect 



on the calculated masses of the corresponding s and d 
states. Both features can be understood by realizing that 
di-electron widths are sensitive to the wave functions at 
the origin whereas spectral masses are more related to 
their long distance behavior. 

B. Bottomonium 

A parallel analysis to the one just carried out for char- 
monium can be done for bottomonium. As ns states up 
to n = 4 have been experimentally identified we choose 
T(3s) and T(4s) to fix the parameters \b and mi,. From 
M [T(3s)] = 10355 MeV and M [T(4s)] = 10579 MeV 
we find Ah = 102.6 MeV.fm and = 4795.5 MeV. The 
results for the spectrum are shown in Table V and as- 
signed to data. For comparison results from the QCDSA 
are also listed. For the sake of completeness it should 
be pointed out that the results for ns states with an 
"intermediate" model based on the asymptotically con- 
stant screened potential of Eq. (|4]) ^1^ lye in between ours 
and the QCDSA ones. On the other hand quark poten- 
tial models not incorporating screening [isj predict much 
larger energy splittings for high n^. 

Again the Ip, 2p and Id states are not well described. 
Now for Ip (2p) with < >i/2= o.4 fm (0.7 fm), the 
discrepancy goes from 60 MeV (30 MeV) for Xb2 to 110 
MeV (70 MeV) for Xfco- For Id, with < >i/2= o.6 
fm, the calculated mass differ about 50 MeV from the 
only known experimental candidate. Unfortunately we 
have not at disposal data for higher radial p or d excita- 
tions to fix a value for the rms radius beyond which spin 
dependent contributions are negligible. If we assume a 
correct prediction for the 2>p state this radius would be 
of ~ 1 fm. 

A very good correspondence between calculated and 
experimental masses (difference of 15 MeV at most) is 
found for Is, 2s, 6s and 7s if the T(10860) is assigned 
to T(6s) (not to T(5s) as usually done) and T(11020) to 
T(7s). Notice that the recent measurements by BaBar 
give 10876 MeV and 10996 MeV for the masses of these 
two resonances. Moreover the T(11020) appears in [2l| 
as a peak between 10990 MeV and 11060 MeV what is 
compatible with being the overlap of our 7s and 8s states. 

It should be emphasized that the assignment of 
T(10860) to T(6s) implies the existence of a T(5s) reso- 
nance with a mass 

M [T(5s)] - 10748 ± 15 MeV (15) 

what can be considered as a main prediction (the quoted 
error of 15 MeV has been estimated from Table V) and at 
the same time as a stringent test of our potential model. 
The presence of this resonance might have some rela- 
tion with the experimental shoulder present on the tail 
of T(4s) with a mass of 10684 ± 10 ± 8 MeV and a width 
of 131 ± 27 ± 23 MeV in reference [221 and a mass be- 
tween 10670 MeV and 10730 MeV in [2i| (see Table I of 
this reference). In the recent study by BaBar [13] there 
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{ML,nr)pDG 


(Ml,„Jsa 




fm 


MeV 


MeV 


MeV 


Is 


0.2 


9458 


9460.30 ± 0.26 


9453 


2s 


0.5 


10037 


10023.26 ± 0.31 


10010 


Id 


0.6 


10218 


10161.1 ± 1.7 


10144 


3s 


0.8 


IO355I" 


10355.2 ± 0.5 


10356 


2d 


0.8 


10471 




10446 


4s 


1.1 


IO579I" 


10579.4 ± 1.2 


10630 


5s 


1.4 


10748 




10862 


6s 


1.7 


10880 


10865 ± 8 
10876 ± 2^" 


11067 


7s 


2.0 


10986 


10996 ± 2^" 
11019 ±8 


11240 


8s 


2.4 


11073 






9s 


2.7 


11144 






10s 


3.1 


11205 






lis 


3.6 


11256 






Ip 


0.4 


9970 


Xw (9912.21 ±0.57) 
Xm (9892.78 ±0.57) 
Xi,o(9859.44 ±0.57) 


9884 


2p 


0.7 


10300 


X62 (10268.65 ±0.72) 
X61 (10255.46 ±0.72) 
Xbo(10232.5 ±0.9) 


10256 


3p 


1.0 


10535 




10541 



TABLE V: Calculated 66 masses and rms radii from V{r) 
with Xt, = 102.6 MeV.fm and rub = 4795.5 MeV. Notation as 
in Table IV. The superindex Ba indicates now recent BaBar 
data 20]. Masses calculated in the QCD string approach are 
taken from [13|. 

appears a small bump around 10700 MeV not identified 
as a resonance (see Fig.l of this reference) that might 
have to do with the predicted state. It should be added 
that the presence of the close BgBs threshold at 10732 
MeV may complicate the experimental extraction of this 
resonance, if existing. 

An additional argument in favour of this resonance can 
be elaborated from the comparison of the experimen- 
tal energy differences between contiguous s excitations 
in charmonium and bottomonium as done in Table VI. 

The assumption that the 4s — 3s mass differences in 
bottomonium (224 MeV) and charmonium have close val- 
ues as it is the case for 3s — 2s and 2s — Is requires for 
cc a 4s resonance around 4260 MeV such as our model 
predicts. Then, assuming that X(4260) is the 4s sate, 
the extension of the argument to the 5s — 4s mass differ- 
ences (158 MeV in charmonium) implies the existence of 
T(5s) about 10740 MeV. Alternatively, as it is the case 
in the QCDSA, the X(4260) could be not a cc state and 
the T(10748) not exist but in such a case the energy dif- 
ference pattern in charmonium and bottomonium would 
be very different (the 4s — 3s mass difference in charmo- 
nium would be 382 MeV against 224 MeV in bottomo- 
nium). It should be emphasized that this discrepancy in 



(n,. + 1) (Mo,„^+i - Mo,„,.) (Mo,n^+i - Mp.n,,) 
cc 66 

2 589 563 

3 353 332 

4 224* 224 

5 158* 
6 

TABLE VI; Experimental mass differences (in MeV) between 
{ur + l)s and n^s states in charmonium and bottomonium. 
The superindex * indicates that the corresponding difference 
have been calculated assuming that X(4260) and ■0(4415) are 
the 4s and 5s states of cc. 

the interpretation of the experimental states is directly 
related to the different manner at which string breaking 
is implemented in both models. Then the experimental 
confirmation (refutation) of our results would serve to 
establish the coulombic (non-coulombic) character of the 
asymptotic quark-antiquark potential. 

Concerning other ns states with n > 8 the small sep- 
aration in energy between neighbors suggest important 
overlaps among them and difficulties for a separated iden- 
tification. This may explain the non-identification of any 
clear signal for a resonance in the region 11000 — 11200 
MeV recently explored by BaBar [20] . Let us realize that 
the limiting mass of the spectrum is still quite far above 

{ML^,n^t)bb ^ nib + m-^ + uv = 11720 MeV (16) 

It is also worthwhile to emphasize the non-relativistic 
character of the fitted static spectrum in bottomonium 
since (pb/mi,)'^ ~ 0.1 — 0.06. Then we can tentatively 
identify Afc = 4(q;s)^/3 being the quark-quark-gluon 
coupling at the bottomonium scale. This gives {as)i^ = 
0.39 in agreement with the value derived from QCD in 
bottomonium for the Ip and 2p states . 

As we are dealing with a non-relativistic system 
we expect that the wave functions obtained from the 
Schrodinger equation may accurately give account of 
other observables. In particular s— wave splittings and 
leptonic (di-electron) widths depend directly on the val- 
ues of the wave functions at the origin. Thus in first 
order perturbation theory the splitting energy between 
the triplet T(ns) and the singlet rjhins) spin states is 
given by 

M [T(ns)] - M Mns)] - ^ (a,), ^ |i?„..(0)l' (17) 

where i?n^s(0) stands for the radial wave function at the 
origin for T(ns). The resulting splitting for n,. = 1 is 

M [T(ls)] - M [7/h(ls)] = 173 MeV (18) 

in accord with the experimental value 

M [T(ls)],, - M [r7fc(ls)]e. = 160 ± 40 MeV (19) 
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For rir = 2 the predicted value is 

M [T(2s)] - M [r]b{2s)] = 70 MeV (20) 

for jirS states 



Regarding the leptonic widths 
they can be evaluated as Q 



1 - 



A(7i,.s) 



(21) 

The terms with (— 16 (as);, /Stt) and (A(rirs)) give ac- 
count of the leading order radiative and higher order ra- 
diative -|- relativistic corrections to 



.(0) 



_ {rirS) 



AeW 

^^b"^ I p 



.(o)r 



(22) 



where ~ ~l/3 is the quark electric charge, a — 
1/137.036 the fine structure constant and M^^s the mass 
of the rirS state for which we shall use the experimen- 
tal value. The calculated leptonic widths are shown in 
Table VIF Although the correction A depends on the 
particular rirS state we shall consider it, for the sake of 
simplicity, as an effective constant. We fix its value from 
Fe+e- (T(10580)) since it corresponds to the highest exci- 
tation with well identified quantum numbers (4s) and we 
expect the non-relativistic and static approaches to be 
more accurate for it than for lower excited states. Then 
we get A = 0.22, one third of the value of the first order 
radiative correction 16 (as)^ /37r = 0.66. 









Is 


1.7 


1.340 ± 0.018 


2s 


0.61 


0.612 ±0.011 


3s 


0.39 


0.443 ± 0.008 


4s 


0.27''^ 


0.272 ± 0.029 


5s 


0.21 


0.20 ±0.05 ±0.10* 


6s 


0.16 


0.22 ±0.05 ±0.07* 


7s 


0.13 








0.095 ±0.030 ±0.035* 


8s 


0.11 





TABLE VII: Leptonic widths r^+^- (in keV) for bb. The 
superindex f indicates the value used to fix A. Data from 
01 except for rirS with rir > S taken from [S^l and indicated 
by a superindex *. The experimental number between the 7s 
and 8s states indicates that the resonance T(11020) in [22|, to 
which this number is assigned, could be a result of the overlap 
of the 7s and 8s states. 

It should be pointed out that the measured 
Fg+g- (10860) in [3| might be contaminated by the hid- 
den T(10748). Instead for > 5 data from 22] where 
a resonance about 1700 MeV is taken into account are 
used. 

A look at the Table makes clear that a good agreement 
with data (10% of difference at most) may be achieved ex- 
cept for Fg+g- (Is). This might have to do either with the 



A dependence on the n^s state or with a deficient descrip- 
tion of the wave function at the origin for T(ls), the more 
relativistic state for bottomonium with the more impor- 
tant spin-spin correction. Indeed a 13% of decrease in the 
value of |i?is(0)| would fit the central experimental value 
of Fg+g- (Is) (notice that the Is spin splitting would be 
130 MeV, still within the experimental uncertainty). It 
should be also kept in mind that a systematic deviation 
of the values of the wave functions at the origin might be 
hidden through the effective value of A. Actually the val- 
ues we get for i?is(0) and i?2s(0) are significantly bigger 
than the ones obtained in the QCDSA. 

Di-electron widths can be also calculated for nd states 
but no data are available. Therefore we will only men- 
tion that the calculated values in our model from (see for 
instance [Tsj ) 



(nd) 



KM 



(23) 



where R^^iO) stands for the second derivative of the ra- 
dial wave function at the origin, are four order of magni- 
tudes smaller than for n^s states. 

For the sake of completeness El decay widths are also 
evaluated. By using a single quark operator approxima- 



tion the width can be written as 25 1 



-rEi 



27 



elakff{2Jf + l)D}f 



(24) 



where kif is the photon energy or momentum, Jj the 
total angular momentum of the final meson and D^f the 
transition matrix element 



Dif = / drui{r)- 



■Jo 1 -TT- ] ~Jl\ -TT ) 

(25) 

being Uij{r) the reduced radial wave functions of the 
initial and final mesons and joi ji spherical Bessel func- 
tions. The results obtained for T(2s) — > 7Xhj(lp) and 
T(3s) — > 7Xfc,/(2p) are compiled in Table VIII. It has 

bb 



Transition 


Tei 




T(2s) 




1.7 


1.2 ±0.2 


T(2s) 




2.5 


2.2 ±0.2 


T(2s) 




2.6 


2.3 ±0.3 


T(3s) 


-*7X6o(2P) 


1.9 


1.2 ±0.2 


T(3s) 


--7X6i(2P) 


3.2 


2.6 ±0.5 


T(3s) 


^7Xfc2(2P) 


3.5 


2.7 ±0.6 



TABLE VIIL El decay widths for bb (in keV) as compared 
to data from Q]. 

to be reminded that the three calculated Xbj{'>T'P) states 
are degenerate in our model. Hence the same wave func- 
tion is employed for all of them. This can be justified 
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by assuming that the experimental masses are explained 
by the effect of spin dependent interactions calculated in 
perturbation theory to the first order. The differences 
in Table Vlll among the three T(2s) 7Xf)j(lp) or the 
three T(3s) jxbji'^p) decays come from the use of the 
non-degenerate experimental XbJ masses to evaluate ki / . 

A clear bias of the results is observed: they are sys- 
tematically higher than data being the discrepancy more 
pronounced for Xboin-p) final states. As Xboinp) states are 
the ones requiring a bigger spin-dependent mass contri- 
bution in our model, the systematics may be suggesting 
the need of implementing Xbjin-p) wave function correc- 
tions. 



V. SUMMARY 

To summarize, a universal form for the quark- 
antiquark static potential, incorporating the screening 
of the color charges by sea pairs, has been proposed 
within a Non-Relativistic Quark Model framework. This 
potential, with a confining long-distance coulombic be- 
havior, reproduces the highly excited light-quark meson 
spectrum and provides a successful spectral description 
of charmonium and bottomonium suggesting the assign- 
ment of A (4260) to the 4s state of cc and the existence 
of a non-cataloged T( 10748) resonance corresponding to 
the 5 s state of bb. These very distinctive predictions of 
our model come from the way screening have been im- 
plemented in it. Therefore their experimental confirma- 
tion or refutation could allow to establish the coulombic 
or non-coulombic character of the long-distance quark- 
antiquark static potential. 

It should be remarked that the only dependence of the 
potential on the particular meson sector comes from the 



value of an effective "gluonic" strength. As the light- 
quark mesons and to a lesser extent the charmonium are 
clearly relativistic systems one can tentatively think that 
some relevant relativistic corrections could be effectively 
taken into account through the value of this parameter. 
The fact that the "gluonic" coulomb strength gets sys- 
tematically a greater value than the gluonic chromoelec- 
tric strength in QCD seems to point out in this direc- 
tion. For the non-relativistic bottomonium this "gluonic" 
strength can be put in correspondence with the strength 
of the chromoelectric one gluon exchange interaction in 
QCD or, equivalently, with the quark-quark-gluon cou- 
pling as at the corresponding scale. 

In our non-relativistic treatment the quark and anti- 
quark masses are parameters of the model. Their values 
should be added to the binding energies to get the meson 
masses. A peculiarity of our potential is the absence of 
any additive constant to get acceptable values of the con- 
stituent quark masses (in the sense of being able to give 
account of other observables such as hadronic magnetic 
moments) from the fitted meson masses. 

All these features make the effective non-relativistic 
quark model proposed very useful to identify excited 
states from existing experimental candidates and for as- 
signing quantum numbers to them. Furthermore it can 
be used to advance predictions on highly excited states 
in all meson sectors. 
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